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Abstract 

Understanding the pattern of gluon radiation in tt production processes is 
important for making an accurate determination of the top mass from the mea- 
surement of its decay products. In a recent paper we showed that the exact 
matrix element and parton shower (HERWIG) calculations gave very different 
results for the distribution of gluon jets in tt production at the Tevatron pp col- 
lider. By repeating the calculation for the simpler e + e~ — * bbW + W~ g process, 
we reveal even more dramatic differences between the two approaches. We con- 
clude that there are significant differences in gluon radiation between HERWIG 
and the matrix element calculation in regions of phase space where one would 
expect agreement. 



Having established the existence of the top quark through its production in pp 
collisions at the Fermilab Tevatron [^], the CDF and DO collaborations must now 
accurately measure the top mass. To do so requires understanding gluon radiation 
in top events. Because top is observed only through its decays, measuring its mass 
requires reconstructing the momenta of its decay products and comparing measured 
invariant mass distributions with those predicted by theory. And because gluons can 
be radiated in top production and decay, top events often contain extra jets which 
may or may not need to be counted among the decay products. While this cannot be 
decided on an event-by-event basis, in principle the problem of extra jets can be dealt 
with as long as their distributions are properly understood. 

In previous work || || we studied the distributions of extra jets in ti production and 
decay at the Tevatron pp collider. In Q] we performed a complete tree-level calculation 
of the process pp — *> bW + bW~+ jet +X, in which we included all contributions from 
gluons emitted in the production and decay stages as well as all top width effects and 
spin correlations. We compared our exact matrix element results to those obtained 
using the parton-shower Monte Carlo program HERWIG ||, which is widely used 
in the experimental analyses. We found a significant discrepancy between the two 
calculations. Although we were unable to identify the reason for the discrepancy, the 
differences in jet distributions seemed to indicate a relative lack of gluons radiated in 
top decay compared to top production in HERWIG. 

In this paper we pursue further the comparison of matrix element (ME) and parton 
shower (PS) calculations of gluon radiation in top production and decay. In light of the 
subtle issues associated with making ME-PS comparisons, and of the complications 
associated with hadronic top production, it is necessary to investigate the differences 
in a systematic way. Therefore we focus here on the simpler case of e + e~ collisions, 
where many of the complications of pp collisions, such as initial state radiation, are 
absent, and where a cleaner comparison can be performed. Then, armed with our 
e + e~ results, we return to a discussion of the hadronic case. 

Comparing a fixed order calculation with a full parton shower program presents 
several challenges. Most notable is that the ME calculation produces only partons 
directly from the hard process, whereas in the PS calculation, any number of partons 
can be present as a result of showers, and multiple gluon effects can be important. 
In order to develop a meaningful method for comparison, we begin by studying the 
well- understood process e + e~ — > qqg, where q is a light quark. Next, we make the 
same comparison but including a non-zero quark mass (e + e~ — > QQg). From these 
processes a useful method of comparison is found, and the regions of applicability for 
the different calculations are studied. Having understood the simpler processes, we 
then consider e + e~ — > bW + bW~ + g via top pair production. Cuts are kept to the 
minimum necessary to avoid singularities. 

In each case, the final state hadrons in the ME calculation are two quarks and 
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a single extra gluon (we ignore the W's from top decay; they can be considered to 
decay leptonically) . In the PS calculation the final state contains two quarks plus any 
number of showered partons. In order to make the ME-PS comparison, we must use 
some algorithm to combine showered partons in the PS calculation so that we are left 
with at most three jets.[] We use the Durham successive combination algorithm ||, 
which has the advantages |?[] that it reproduces LEP data well and lends itself to direct 
comparison with ME calculations. 

We apply the Durham algorithm as follows. For each pair of partons (or jets) i 
and j, we compute the quantity 

2min(E 2 ,Ej)(l-cos%) 
Vij = ~ > I 1 ) 

where is the energy of parton Qij is the angle between them, and s is the 
process center-of-mass energy. The pair with the smallest y^ are combined into a 
single jet with four- momentum p = pi + pj. The jet replaces the partons in the 
Uij computations and we repeat the process until no more than three jets remain. 
Note that for small angles, the numerator of j/y is simply the fcf. of the less energetic 
parton with respect to the more energetic one. Hence the algorithm causes the partons 
with the smallest relative transverse momentum to be combined. In addition, the 
dimensionless variable y is somewhat like an angular variable, so that this method is 
analogous to the angular cone algorithm used in our comparisons for pp |4| . 

After using the jet algorithm, we have three-' jet' events from both PS and ME 
calculations, and they can now be compared. To do so, we make further use of the 
Durham variable: for each event, we plot the distribution as a function of the smallest 
value of y (= y m in) for all pairings of the remaining three jets. Because large values 
of 2/min correspond to large relative k T , i.e., large angles and energies, the ME and 
PS results should agree well in this region. Small values of y min bring us into the 
soft /collinear regime, where the ME result diverges as y m i n — > 0. In this region multiple 
gluon effects become important, and we expect the PS result, which takes these into 
account via resummation, to remain finite. 

This behavior is illustrated in Fig. 1(a), which shows the y m { n distribution for 
e + e~ — > qqg (massless quarks) with a center-of-mass energy of 100 GeV. Since we are 
primarily interested in shapes of distributions, the normalization of the ME calculation 
is chosen (via the ot s scale) so that the distributions for e + e~ — > qqg agree at large 
values of y m i n . We see, then, that ME and PS agree quite well for large y m i n , not 
only in normalization, which we have fixed, but in shape, which we are not free to 
manipulate. As y min decreases, the ME distribution rises, diverging as a s ln 2 (l/y min ), 
but we see the multigluon effects begin to show in the PS curve, which turns over and 
remains finite. 

1 The PS results were obtained using HERWIG || (v5.8), with hadronization and b quark decays 
turned off. 
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The effect of giving the final quark a mass is shown in Fig. 1(b) for the process 
e + e~ — > bbg. Note that the b mass has weakened the ME divergence, and the distri- 
bution for all 7/min is less steep than in the massless case. The PS calculation again 
reproduces the ME result (including the rounding due to the mass) quite well in the 
large y min region, and we see the two results begin to disagree as multigluon effects 
come in for smaller y m i n . 

Satisfied that we understand the above simple cases, we turn to top production. 
The ME and PS calculations performed here for e + e~ production of top are similar 
to those performed in ||. We take m t = 175 GeV. To obtain the ME results we 
have performed an exact parton-level calculation of e + e~ — > bW + bW~ + g, including 
top width effects, using helicity amplitudes generated by the MadGraph || package. 
The PS results were obtained using HERWIG || (v5.8), with hadronization and b 
quark decays switched off. As in the case of hadronic top production, we do not 
consider radiation off the W decay products, and the W's can be considered to decay 
leptonically. 

We proceed with the jet algorithm as above, and compare the ME and PS y min 
distributions for center-of-mass energy 420 GeV in Fig. 2.Q The difference is quite 
dramatic. The ME curve displays further rounding as one might expect from a higher 
mass quark or from the lack of correlation between the final state b quarks and gluons 
radiated from the t or i. The behavior of the PS curve is unexpected. We see the 
usual turn-over at small values of y m i n , but the large y m i n behavior is surprising. Not 
only does PS not reproduce ME; its shape is so different that there is no shift in 
normalization that would cause the two to agree. 

Given that the matrix element calculation is of fixed order in perturbation theory, 
one expects that its region of validity would start at large values of y m i n and extend to 
the left to the region where we expect multiple gluon emissions to become important. 
Since the parton shower has these multiple emissions, but uses an angular ordering 
algorithm, we expect its region of validity to be from small values of y m { n and extend 
to the right. The fact that there is almost no region of overlap is very disturbing. 

We have been unable to pin down the reason for the above discrepancy. Although 
it is in principle difficult to compare fixed order (ME) with approximate all-orders 
(PS) calculations, the fact that the comparison works well for light quarks suggests 
that this is not the explanation for the difference. Furthermore, there is no reason 
why a matrix element calculation which works well for light quarks should fail for 
heavy quarks, provided of course that appropriate cuts are used to avoid the soft and 
collinear regions. Following the results of Ref. Q, it has been suggested that the 
observed differences between the PS and ME calculations were due to a choice of cuts in 
the latter that were too sensitive to the infra-red region, thus implying that important 

2 This energy is chosen to approximate the average subprocess center-of-mass energy for ti pro- 
duction in pp collisions at the Tevatron. 
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higher-order multi-gluon contributions were missing from the ME calculation. However 
(see below) increasing the jet E™ in cut in the calculation of Ref. |4j], and thus moving 
to a phase-space region where the ME approach should work even better, does not 
resolve the problem. 

We are therefore forced to conclude that the implementation of very heavy quark 
production and/or decay in HERWIG may not be correct. One can perhaps shed 
further light on this by comparing the two calculations for stable top production, 
thus testing the 'production' part of the jet cross section in each case. The result 
of such a comparison is shown in Fig. 3, where the distribution in y min is shown for 
stable m t = 175 GeV tig production in e + e" collisions at \/s = 420 GeV. In contrast 
to Fig. 1 (essentially the same distribution but for light quarks), there is a marked 
difference between the two calculations in the medium-small y m \ n region. Note that 
in this case there is no collinear singularity and so the divergence as y m \ n — > is 
weaker (~ a s \ny min ) in the ME distribution. It is difficult to see how higher-order 
multi-gluon contributions could produce such a pronounced increase at y m m ~ 10~ 4 
in the PS calculation. By comparing Figs. 2 and 3, we see that the difference in the 
'production' part can also be seen in the full 'production + decay' distribution. 

Finally, we return to the case of hadronic top production. We recall that the 
discrepancy we found in M is similar to that observed in e + e~ collisions (Fig. 2). Es- 
sentially it amounts to a relative deficit of decay stage radiation in the PS calculation. 
This effect is both enhanced and clarified in the study in e + e~ because of the lack of 
initial state radiation, and the cleaner environment. As a final confirmation of this 
fact, we rep lot in Fig. 4 the distributions in the phase-space separation variable Aitk- 
in pp collisions at yfs = 1.8 TeV from the ME and PS (HERWIG) calculations, with 
a very large cut on the transverse energy of the jets. In this region, the tree-level 
ME calculation should certainly be a good approximation [S]. Again we observe a 
large discrepancy between HERWIG and the matrix element calculation. It is difficult 
to estimate the effect of this difference on top quark analyses, and in particular on 
the measurement of the top mass, but clearly every effort should be made to ensure 
that the event simulators used to study the top quark do indeed constitute a good 
approximation to the underlying physics. 
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Figure Captions 

[1] (a) e + e~ -> qqg at = 100 GeV: PS (HERWIG) versus ME as a function of 
7/ m i n . A cut of y ni \ n > 2 x 10~ 3 has been placed on the matrix element calculation, 
(b) e+e- -> big at ^/s = 100 GeV: PS (HERWIG) versus ME as a function of 
2/min- A cut of 2/ m i n > 2 x 10~ 4 has been placed on the matrix element calculation. 

[2] e+e" -> bW + bW~g at y/s = 420 GeV with m t = 175 GeV: PS (HERWIG) 
versus ME as a function of y m i n . 

[3] e + e~ — > tig at ^/s = 420 GeV with m t = 175 GeV and no top decay: PS 
(HERWIG) versus ME as a function of y m in- 

[4] pp -> bW + bW-g + X at ^/s = 1.8 TeV with m t = 175 GeV: PS (HERWIG) 
versus ME, distribution in AR bj for Ef n > 20 GeV, AR bj > 0.4. 
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